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We analyse the discovery prospects of the charged Higgs boson, H± , via its decay in the W±h 
channel in the minimal supersymmetric Standard Model (MSSM) as well as several 2-Higgs Doublet 
Models (2HDMs). h , the lightest scalar Higgs boson in these models, is identified with the recently 
discovered ~ 125GeV state, H 0 b s , at the Large Hadron Collider (LHC). We find that, while it 
provides an important input in the kinematic selection of signal events, the measured H 0 b s mass 
renders this channel inaccessible in the MSSM. In the 2HDMs though, through a dedicated signal- 
to-background analysis for the pp —y t(b)H~ —>• T^vijjbbbtf))-^ h.c process, we establish that some 
parameter space regions will be testable at the LHC with yfs — 14 TeV. 


I. INTRODUCTION 

The MSSM and the 2HDMs contain two scalar electroweak (EW) doublets, 4>i and 4>2, which results in a total 
of five physical states, two scalars, h , H , one pseudoscalar, A, and a charged pair, H ± , upon EW symmetry 
breaking. The mass of the H ± is a free parameter in these models. When the H ± is heavier than the top quark, 
its main production process at the LHC is pp —> tH~b (by which we actually imply pp -> tH~b + pp — tH + b] 
here and onwards, we do not distinguish between fermions and anti-fermions when their identity is unspecified 
and/or can be inferred from the context), and it decays dominantly into tb. However, this decay channel is very 
difficult to probe when the parameter tan /3 (the ratio of the vacuum expectation values of the two Higgs fields) 
has a value less than ~ 3, owing to the large reducible and irreducible backgrounds. For such tan/3, earlier 
studies pQ have shown that the decay channel > W ± h could prove crucial in the MSSM. 

Since the discovery of iL 0 bs at the LHC 00 , its mass measurement, ra# obs ^125 GeV, provides an additional 
input for collider analyses as well as an important constraint on physics beyond the Standard Model (SM). We 
therefore revisit the -D W^h decay channel, assuming h to be consistent with the H 0 b s , in the MSSM, the 
2HDM of Types I and II (inferences drawn for which can largely be extended, respectively, to Types III and IV 
also), and the aligned 2HDM (A2HDM). In all these models, instead of the h , the H can also play the role of 
the iL 0 bs, but we do not consider this scenario here. We estimate the signal strength of the considered process in 
the considered models’ parameter space regions satisfying some important experimental constraints. We then 
compare these estimates to the expected sensitivity obtained from a signal-to-background analysis for the LHC 
with y/s = 14 TeV. 


II. MODELS AND PARAMETERS 


As noted above, the production process of interest here is pp tH~b. When calculating the cross section 
for this process, the two subprocesses, gg tH~b and gb —)> tH ~, cannot be added naively, as that results 
in a double counting between the two contributions. This is because the gg amplitude can be seen as a tree- 
level contribution to the next-to-leading order amplitude that includes a virtual 6-quark, with the gb process 
making the leading order amplitude. The square of the amplitude for the gb — tH~ process, with spin/colour 
summed/averaged, is given by [4] 


\M\ 2 


9_Ih± 2 0»-rag±) 2 

2m^, 4 N c tb s(m 2 — t) 


1 + 2 


m 2 H± - 
u — m 2 H± 



t-m\ 


TO g± ) 

u - m 2 H± ) 


( 1 ) 


where g s is the SU(3)c and #2 the SU(2)l gauge couplings, Nc = 3 is the number of colours and Vtb is the 
relevant CKM matrix element. The model-dependence of this process is embedded in the coupling g q H ±, along 
with the predicted H ± mass. The H ± —» W ± h decay channel, which is proportional to the coupling 

9hH+w- = y (cos/ 35'12 - sin /3S n ), (2) 

where Sn and Si 2 are the elements of the scalar Higgs mixing matrix, has a strong dependence on tan/3, 
similarly to the production process, as we will see below. 
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In the MSSM the squared H^-quark coupling noted above is given by g qH ± = m 2 tan 2 /3+m 2 cot 2 /3, where 
and m t are the masses of the bottom and top quarks, respectively. In addition to the most relevant parameters 
here, tan f3 and m#±, some other free parameters of the MSSM are also crucial, particularly for obtaining the 
desired mass of h. These include the Higgs-higgsino mass parameter, /q the soft masses and trilinear couplings 
of the sfermions and the soft gaugino masses. (These supersymmetric (SUSY) parameters are taken into account 
here under certain unification assumptions; see [5] for further details about them and their scanned ranges.) 

In the 2HDMs, the physical masses of the five Higgs bosons as well as the fermion Yukawa couplings are 
all free parameters. However, in order to avoid large tree-level flavour-changing neutral currents (FCNCs), a 
^-symmetry is generally introduced, so that each type of fermion only couples to a single doublet [HE!. In 
2HDM Type I (2HDM-I), the fermions only couple to <f >2 and in Type II (2HDM-II), the down-type quarks 
and leptons couple to <I>i and up-type quarks to 4>2- The Yukawa couplings are then determined entirely by 
the parameter tan/3. In the CP-conserving limit assumed here, tan/3, | sin(/3 — a)| and m 2 2 are then the only 
additional free parameters. Another mechanism for controlling FCNCs is to require that the two Higgs doublets 
have Yukawa matrices which are proportional to one another, or aligned. Here we consider the case where 
all fermions couple to both <I>i and $2 with aligned couplings, known as the A2HDM. In this model, | shm|, 
A2, A3, | A7 1 and \/3 u,D,L \ are the only free parameters besides the Higgs boson masses. (See [Hjfor a complete 
description of these parameters.) The coupling g qH ± in each of these 2HDMs is given in table hi 


Table I: The expressions for g qH ± in the 2HDMs considered in our analysis. 



2HDM-I 

2HDM-II 

A2HDM 

9qH± 

ml cot 2 /3 + mi cot 2 0 

tan 2 /3 + m 2 cot 2 0 

ml tan 2 + m 2 tan 2 /3 U 


We performed scans of the above mentioned parameter space of each model to search for their regions with 
a strong signal. The particle mass spectra and decay branching ratios (BRs) in the MSSM and the 2HDMs 
were generated using the programs SUSY-HIT-vl.3 [9] and 2HDMC uni, respectively. In these scans we 
required 123 < mh < 127 GeV and 200 < m H ± < 500 GeV, and additionally in the case of the 2HDMs, 
135 < run < 500GeV and tua = m H ±. We allowed a slightly extended range, 1.5 < tan/3 < 6. For the MSSM 
and AHDM, the 5-physics observables BR (B X s ^), BR (B u — » tv), and BR (B s — >> /i + /i _ ) were calculated 
using SuperIso-v3.4, and were required to lie within the 95% confidence level limits suggested in the manual 
of the program DU- For the ^-symmetric 2HDMs, the relevant input parameters were simply chosen so as to 
satisfy these constraints, following [12]. Finally, each scanned point was tested with HiggsBounds-v4.1.3 [13] to 
make sure that the Higgs states other than h satisfy the LEP, Tevatron, and LHC constraints. 


III. SIGNAL-TO-BACKGROUND ANALYSIS 

For our detector analysis, in order to directly reconstruct the H 0 b s (implying here a 125 GeV Higgs state) 
we consider only the iL 0 bs —^ bb decay, since both the 5-quarks are observable and also because in the models 
studied here the h decays dominantly in this channel. The complete process analysed is then pp —)> {b)tH ± 
(b)bW T W ± H ohs —)> ( b)bbbjjlvi , where one of the W-bosons (from either or t-quark) decays leptonically 
and the other hadronically. The presence of a single lepton allows us to avoid multi-jet backgrounds, while 
requiring one hadronic W avoids additional unseen neutrinos. The main background for this process is ttb(b ), 
where either an additional 5-tagged jet combines with a 5-jet from a top decay or an additional 55 pair mimics 
an Hobs —>• 55 decay. 

In order to estimate the sensitivity obtainable at the 14TeV LHC, we generated the t(b)H~ signal using 
Pythia 6.4.28 [14] with the MATCHIG [15] add-on to avoid double counting among the gg —>> tbH~ and 
gb —>• tH~ processes, and all t(b)WX,X —)> bb backgrounds with MadGraph5 [16]. The parton showering and 
hadronisation for both the signal and the background was performed using Pythia 8 m and the subsequent 
detector simulation with DELPHES 3 [18], using experimental parameters based on the ATLAS experiment. 
The 5-tagging efficiency chosen was e ?7 tanh(0.03 pr — 0.4), with = 0.7 for central (\rj\ < 1.2, with \rj\ being 
the rapidity) and = 0.6 for forward (1.2 < \rj\ < 2.5) jets, and the transverse momentum, pt, in GeV. 

To reconstruct our signal events and separate them from the background, we then proceeded as follows: 

1. Only events with > 3 5-tagged jets, > 2 light jets, one lepton (e or fi) and missing energy (J^ T ) > 20 GeV 
were selected. For all objects we required pt > 20 GeV, |^| < 2.5 and A R, their separation from other 
objects, > 0.4. 

2. First the pair of light jets with an invariant mass, nijj , closest to the W boson mass, mw , was chosen. 
An event was rejected if no pair satisfied \rrijj — mw | < 30 GeV. 
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3. The observed lepton was used to find the longitudinal component of the neutrino momentum, p u , z , by 
imposing the mass constraint mf> v — mw , thus attributing all $ T to a neutrino from a W decay. Such 
solutions have a twofold ambiguity as a result of the quadratic nature of the constraint. Both solutions 
were kept when they were real. For complex solutions, a single real p u ^ z was retained, discarding the 
imaginary component. 

4. The background can mimic the signal by combining a 6-jet from a top decay with an additional 6-tagged 
jet to reconstruct the h. Thus as a ‘top veto’ was applied at this stage, but only in the large m H ± region. 
This implies that an event was rejected if two t-quarks could be reconstructed from reconstructed W ’s 
and any unassigned jets, with both satisfying \mwj — m t \ < 20 GeV. The jets used may or may not be 
6-tagged. For m H ± not too far above the H ± —>> W ± H 0 b s threshold, one of the resulting 6-jets combines 
with the W± to give an invariant mass m^w ~ m t in a large fraction of the available phase space. Such 
signal events are cut if the top veto is applied at this point, negating the benefits of the background 
reduction. Therefore, for m H ± < 350 GeV, the top veto was postponed until after the next step. 

5. The pair of 6-tagged jets with invariant mass, ra^, closest to ra# obs ^ 125 GeV was used to reconstruct 
the i^obs- The event was rejected if no pair satisfied | — m# obs | < 15 GeV. 

6. In the low m#± region, the top veto was applied at this stage in the same way as above, but excluding 
the 6-jets used in H 0 b s reconstruction. 

7. From the reconstructed IT’s and remaining 6-tagged jets, the best top quark candidate was identified 
as the one for which the Wb combination had an invariant mass, mwb , closest to m t . If the selected 
combination included one leptonic W solution, the other was discarded. The event was rejected if there 
was no good candidate with | mwb — m t \ < 30 GeV. 

8. The reconstructed H 0 b s was combined with the remaining W to yield the discriminating variable mwH ohs - 
If there were two leptonic IT’s remaining, both values of mwH ohs were retained. For each m H ± considered, 
we placed a cut on the range of reconstructed mwH ohs which maximises the statistical significance, S/y/~B. 


IV. RESULTS AND DISCUSSION 

In the MSSM, for the selected m H ± range, mh lying above 123 GeV can only be reached (for 0(1 TeV) scanned 
values of the soft SUSY parameters) near the allowed upper limit of tan / 3 , as shown in the left panel of Fig. [l] The 
right panel shows that such values of tan (3 and always yield a very poor a(pp —)► ti7 ± ) x BR(H ± W ± h). 
Multiplying it with BR (h —» 66) gives the total signal strength for a given point, which is further reduced, as is 
evident from the heat map in the panel. We therefore do not test this model against the collider analysis. 



Figure 1: m h as a function of m H ± and tan f) in the MSSM (left) and cr(pp tH ± ) x BR (H ± —>> W±h) as a 

function of m H ±, with the heat map showing the BR (h 66) (right). 

For the 2HDMs, Fig. [2] shows the results of the parameter scans along with the sensitivity expected from the 
collider analysis. For the ^-symmetric cases, we find a large number of points which are potentially discoverable 
at a high-luminosity (3000 fb -1 ) LHC, as seen in the left panel. The A2HDM shows even stronger signals, well 
within reach of even the standard luminosity (300 fb -1 ) LHC. In the right panel we consider the effect of 
imposing the signal strength, p x , where p x — a(pp H 0 b s — > X)/a(pp — >• Hsm —> X), with Hsm being a 
125 GeV SM Higgs boson, measurements on h for these points. For this purpose, the theoretical counterparts of 
fi x (with i7 0 bs replaced by h) were determined with HiggsSignals-vl.20 [T9\ for X = 77, ZZ and compared with 
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the CMS measurements, p 1 11 = 1.13 ±0.24, p zz = 1.0 ±0.29 [20] . We note that, owing to significant deviations 
of h from iLsM-like properties for most of the good points from the scans, those surviving after imposing the 
p x constraints have a significantly lower detection potential, particularly in the case of ^-symmetric 2HDMs. 
In the A2HDM, some points still remain testable even with 300 fb -1 luminosity. We therefore conclude that 
the H ± —>• W ± h channel can be a useful probe of in the 2HDMs with high luminosities at the LHC. 
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Figure 2: Signal strength for the points satisfying only the 6-physics constraints (left) and additionally the 
CMS constraints on p 11 and p zz (right) from the scans for the 2HDMs. Also shown are the contours for 
different statistical sensitivities expected for selected integrated luminosities at the LHC with y/s = 14 TeY. 
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